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Abstract 

A theoretical analysis of the partial inclusive B —> X^y decay rate with a cut Ey > Eq on photon 
energy must deal with short-distance contributions associated with three different mass scales: 
the hard scale mb, an intermediate scale and a soft scale A, where A = mt- 2Eo ~ 1 GeV 

for £0 w 1.8 GeV. The cut-dependent effects are described in terms of two perturbative objects 
called the jet function and the soft function, which for a next-to-next-to-leading order analysis 
of the decay rate are required with two-loop accuracy. The two-loop calculation of the soft 
function is presented here, while that of the jet function will be described in a subsequent paper. 
As a by-product, we rederive the two-loop anomalous-dimension kernel of the B-meson shape 
function. 



1 Introduction 


Weak-decay processes involving flavor-ehanging neutral eurrents are sensitive to the effects of new 
physies, beeause the deeay amplitudes are loop-suppressed in the Standard Model. In this eontext, 
the deeay B Xsj plays an espeeially prominent role, sinee its rate is being measured inereasingly 
well by the 5-faetories. The eurrent experimental preeision already matehes the theoretieal aeeuraey 
of the next-to-leading logarithmie predietion. This has triggered an effort to push the preeision of 
the theoretieal ealeulation of the deeay rate in the Standard Model to the next level of aeeuraey. Due 
to the presenee of several different seales in the deeay proeess, this ealeulation involves a number of 
different elements. 

Several of the steps required to aehieve next-to-next-to-leading logarithmie (NNLO) aeeuraey 
have already been taken. The matehing of the Standard Model onto an effective weak Hamilto¬ 
nian has been eompleted by performing a three-loop matehing ealeulation onto the eleetro- and 
ehromo-magnetie dipole operators ||T||. The effeetive weak Hamiltonian allows one to resum large 
perturbative logarithms of the form Us \n(Mwlmb). To this end, the three-loop anomalous-dimension 
matriees for the four-quark operators [2] and for the mixing of the dipole operators among eaeh 
other have been ealeulated. The evaluation of the four-loop anomalous dimension of the mixing 
of the eurrent-eurrent operators into the dipole operators is in progress and is now the only missing 
element to obtain the Wilson eoefheients in the effeetive weak Hamiltonian with NNLO preeision. 
The most difheult part of the ealeulation is the evaluation of the matrix elements of the eorrespond- 
ing operators, in partieular the ones involving penguin eontraetions of eurrent-eurrent operators with 
eharm-quark loops. So far, the eomplete two-loop matrix element is known only for the eleetro- 
magnetie dipole operator Qjy 0|. For other operators, only the parts proportional to jSoa? (more 
preeisely, the terms proportional to rifa^) have been obtained Q. For the ease of Qjy, not only the 
total deeay rate but also the photon-energy speetrum has been evaluated at two-loop order For 
the remaining operators in the effective weak Hamiltonian, the photon spectrum is known to order 

ySoa? nil. 

An important motivation for undertaking a NNLO (i.e., order a^) evaluation of the B X^y 
deeay rate is the faet that sueh a ealeulation would reduee the strong dependenee on the renormal¬ 
ization seheme adopted for the eharm-quark mass, whieh is observed at NLO |IH1E|. However, the 
eharm mass is not the only low seale in the deeay proeess. Another set of enhaneed eorreetions 
arises beeause it is experimentally neeessary to put a out Ey > Eq on the photon energy (defined 
in the 5-meson rest frame). The relevant seale is A = nib - '^Eq. For Eq « 1.8 GeV, the currently 
lowest value of the out aehieved by the Belle experiment ifTTll . the seale A « 1 GeV is barely in the 
perturbative domain. For even higher values of Eq, the effeots assoeiated with the seale A eannot 
be calculated reliably in perturbation theory, in whieh ease they are relegated into a nonperturbative 
shape-function IfTTlIT^II . 

As long as the out energy is ohosen sufhoiently low, sueh that » A » Aqcd, the partial 
inolusive B X^y deeay rate oan be ealeulated perturbatively using a multi-soale operator-produot 
expansion. It oonsists of a simultaneous expansion in powers of A/m^, and Aqcd/A, oombined with 
a systematio resummation of logarithms of ratios of the hard, intermediate, and soft seales ifOl 
(see also d). At leading power in A/m^ and next-to-leading order in the expansion in powers of 
Aqcd/A, it is possible to derive an exaet expression for the partial deeay rate r(A), valid to all orders 
in perturbation theory, in whieh the dependenee on the variable A enters in a transparent way. The 
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Here nib is the Z?-quark pole mass, and nibip) denotes the running mass defined in the MS seheme. 
The only hadronie parameter entering at this order is the quantity related to the fi-quark kinetie 
energy inside the B meson. The ellipses represent subleading eorreetions of order (Aqcd/A)^, whieh 
are unknown. The pole mass and must be eliminated in terms of related parameters defined in 
a physieal subtraetion seheme, sueh as the shape-funetion seheme ifT^ [TTl . The seales ~ nib, 
/Xi ~ ylnibA, and //q ~ A are hard, intermediate, and soft matehing seales. The hard funetion Hy, 
the jet funetion j, and the soft funetion 's eneode the eontributions to the rate assoeiated with these 
seales. Note that all information about the short-distanee quantum fluetuations assoeiated with the 
weak-interaetion vertiees in the effeetive weak Hamiltonian are eontained in Hy. Logarithms of 
ratios of the various seales are resummed into the evolution funetions U\ (evolution from the hard to 
the intermediate seale) and U 2 (evolution from the intermediate to the soft seale), as well as into the 
quantity 

77 = 2 r — rcusp[Q'^ 0 u)], ( 2 ) 

whieh is given in terms of an integral over the universal eusp anomalous dimension of Wilson loops 
with light-like segments IfT^ . The result dT]) is formally independent of the ehoiees of the mateh¬ 
ing seales. In praetiee, a residual seale dependenee remains beeause one is foreed to truneate the 
perturbative expansions of the various objeets in the formula for the deeay rate. Redueing the seale 
uneertainty assoeiated with the lowest short-distanee seale, A « 1 GeV, is the goal of the present 
work. 

The soft funetion Tin dT]) is related to the original fi-meson shape funetion S (co, /x) ifTTI through 
a series of steps. Starting from a perturbative ealeulation of the shape funetion in the parton model 
with on-shell Z7-quark states, we first define 

I S parton(^5/^) • (3) 

0 

For Q » Aqcd, this parton-model expression gives the leading term in a systematie operator-produet 
expansion of the integral over the true shape funetion in?ll . The first power eorreetion is linked to 
the leading term by reparameterization invarianee 111^11^ and gives rise to the term proportional to 
1 x 1 / in dD- While the perturbative expression for 5 parton involves singular distributions ifT^ . the 
funetion s has a double-logarithmie expansion of the form 



= 1+1; (^) ( 4 "' + 4">L + ■ ■ ■ + 41 ,. 

n=l ' ' 

The funetion T is then obtained by the replaeement rule ||T?1 

T(L,yU) = 5(L,yU) 




( 4 ) 

( 5 ) 


2 












where /„(jc) are n-th order polynomials defined as 
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By solving the renormalization-group equation for the soft funetion order by order in perturbation 
theory, the eoeflieients of the logarithmie terms in ( 0 ) ean be obtained from the expansion eoeffi- 
eients of the shape-funetion anomalous dimension and the /?-funetion, together with the eoeflieients 
eoefReients arising in lower orders llT5l . The two-loop ealeulation performed in the present pa¬ 
per gives the eonstant c® and provides a eheek on the two-loop anomalous dimension of the shape 
funetion. We also note that from our result for s{L,fj.) one ean derive the two-loop expression for 
SpaitonioJ, fj) in terms of so-ealled star distributions ifT^ . 

In the next seetion, we diseuss how to perform the two-loop ealeulation of the soft funetion s in 
an efheient way. The ealeulation is simplified by representing the d-funetion operator appearing in 
the shape-funetion as the imaginary part of a light-eone propagator. In this way, we avoid having 
to deal with distribution-valued loop integrals and instead map the ealeulation to the evaluation 
of on-shell two-loop integrals with heavy-quark and light-eone propagators. Using integration-by- 
parts relations among these loop integrals, the entire ealeulation is redueed to the evaluation of four 
master integrals. After presenting the result for the bare soft funetion, we diseuss its renormalization 
in Seetion 0 The relevant anomalous dimension depends both implieitly (through the eoupling 
eonstant) and explieitly (through a star distribution) on the renormalization seale. This explieit 
dependenee gives rise to Sudakov logarithms in the soft funetion. We eonelude after presenting our 
final expression for the renormalized soft funetion in Seetion 0 


2 Two-loop calculation of the soft function 


The definition of the soft funetion 5 in © implies that 
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where d,, are effeetive heavy-quark fields in heavy-quark effeetive theory lOTll . are on-shell b- 
quark states with veloeity v, and n is a light-like 4-veetor satisfying n ■ v = \ (note that = 1 and 
= 0). The normalization of states is sueh that (d,,! d,, d,, Id,,) = 1. 

Working with the above representation of the soft funetion is diflioult due to the presenee of the 6- 
funetion differential operator, the Feynman rules for whieh involve 6 funetions and their derivatives. 
This eomplieation ean be avoided by writing the d-funetion operator as the diseontinuity of a light- 
eone propagator in the baekground of the gluon field. This allows us to represent the soft funetion 
as a eontour integral in the eomplex oj plane: 
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The Feynman rules for the gauge-eovariant propagator involve light-eone propagators of the type 
(to + n ■ py^, whieh are straightforward to deal with using dimensional regularization and standard 
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Figure 1: Two-loop graphs contributing to the soft function. Double lines denote heavy-quark 
propagators, while crosses denote possible insertions of the operator {co + in ■ D + iO)"'. 


loop techniques. Dimensional analysis implies that an n-loop contribution to the matrix element in 
the integrand of the contour integral is proportional to where d = A-2eis the dimension 

of space-time. The relevant contour integration yields 
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The two-loop Feynman diagrams contributing to the matrix element in ([HI) are shown in Fig¬ 
ure [H They are on-shell heavy-quark self-energy diagrams with an operator insertion of the gauge- 
covariant light-cone propagator. Instead of drawing a separate diagram for each insertion, we draw 
the topology for a set of diagrams and indicate with a cross the locations where the operator can be 
inserted. The loop integrals arising in the calculation of the soft function contain heavy-quark as 
well as light-cone propagators. The one-loop master integral is 
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where m = a» -l- /O, and 


h{a,b,c) = 


r(a + b- I)F(2a + b + c - d)T{d-2a - b) 


Y{a)Y{b)Y{c) 

The most general two-loop loop integral we need has the form 

-ai-ai-bi -b2-h-Ci -C2 


( 11 ) 


/ 


d^kd‘^1- 


(P)"* (/2)«2 [(k - /)2]“3 (v . kf' (v • [v(k + (n-k + cof' {n-l + 
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where all denominators have to be supplied with a “-l-iO” prescription. Note that we do not restrict the 
exponents ai,..., C 2 to be positive. Loop integrals with non-trivial numerators are written as linear 
combinations of integrals for which some of the indices take negative values. A third light-cone 
propagator, [n • (k - /) -l- cu] “ k can be eliminated using partial fractioning followed by a shift of the loop 
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momenta. We use integration-by-parts identities to reduee the two-loop integrals to a minimal 
set of master integrals. These linear algebraie identities are derived by observing that integrals over 
total derivatives vanish in dimensional regularization. With I 2 = l 2 {ai,a 2 ,aj,,b\,b 2 ,bi„C\,C 2 ), the 
four relations obtained from inserting eaeh of the differential operators dkk^, dk^^n^, and 
into (d) take the form 


0 = [d - 2 a\ - a-i - b\ - c\ - aj -1- aa aj aj - bj, bj -1- c\ c|] h , 

0 = [aa - ai - a\ a| a^ + a\ a| aj - aa a^ aj - 1 - aa a^ aj 
- Z?! bl b2 - Z^a b2 bj -h ci c[ - ci c[ C2] h , 

0 = [a\ a[ - a\ a| cj - aa a^ cj -1- aa C2 -1- Z^i b| -1- Z73 b^] h , 

0 = [ai a\ bj -I- aa aj bj - aa a^ b^ - 2 b\ b| - 2 Z?a b^ - ci c[] h ■ ( 13 ) 


In these equations, the operator a^ (a“) raises (lowers) the index a„ by one unit. Four additional 
relations (with indiees 1 <-^ 2 interchanged) are obtained taking derivatives with respeet to Z^. Also, 
there are partial-fraetion identities for integrals eontaining three different heavy-quark propagators, 
whieh follow from the relation bj -i- b^ - b^ = 0. Repeated applieation of these identities ean be 
used to set at least one of the bi exponents to zero. 

The reduetion to master integrals ean be performed using eomputer algebra. To do so, one gener¬ 
ates the equations for the index range relevant for a given ealeulation and uses Gaussian elimination 
to express eomplieated integrals in terms of simpler ones |E21- Sinee the number of equations is very 
large, the order in whieh the equations are solved is erueial, and ll^ devised an effieient method to 
perform the reduetion. A fast implementation of this algorithm is available in the form of a program 
that solves the relations (fl^ and generates a database eontaining the result for eaeh integral Il2?l . At 
the end of the proeess, we are left with the master integrals 
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(14) 


The first two are evaluated by using the well-known result for the one-loop self-energy integral to 
perform the first loop integration. The result takes the form (fTTib with non-integer exponents, so that 
the seeond loop integration is trivial. The only master integral that needs speeial eonsideration is the 
third one. Its Feynman parameterization is 

Ma = F(7 - 2d) F(f - I)^ f dx f dy(l- x)^-^(I - y)^-y‘^-\l - xy )^-5 . ( 15 ) 
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We evaluate the parameter integral by expanding 
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n=0 


r(2 - f) 

r(n + l)r(2 -n-f) 




(16) 


and summing up the series after the integration. 

With the integrals at hand, the evaluation of the diagrams is straightforward: first, eaeh diagram 
is written in terms of integrals of the form (fT^ . and then eaeh integral is expressed in terms of the 
master integrals. Adding up the eontributions of all graphs and expanding in e = 2 - we find 
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aeeounts for the renormalization of the bare eoupling eonstant, and 
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is the first eoeffieient of the /?-funetion. Throughout, as = as(jS). Note that the bare soft funetion 
is seale independent, sinee the seale dependenee of asipi) eaneels against the explieit ju dependenee. 
Also, sinee the soft funetion is an on-shell matrix element of a gauge-invariant operator it is gauge 
invariant. The terms of 0{e) and O(e^) in the one-loop result above are needed for the evaluation of 
the eounter-term eontributions, as deseribed in the next seetion. The two-loop ooeffieients are 
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3 Renormalization of the soft function 

As usual, we define an operator renormalization faetor Z via 

S(a>,IU)= I doj'Z(a>,CL)',lu)Sbare(Oj'), (21) 
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where Z absorbs the UV divergenees of the bare soft funetion, sueh that the renormalized soft fune- 
tion is finite in the limit 6^0. Here and below, all integrals over variable co, co' ete. run from 0 to oo. 
The eonvolution integral ean be understood as a generalization of the matrix formula O, = Z, , 
and so the usual relation between the Z faetor and the anomalous dimension holds. It follows that 

y(m,= - I doj ---Z (m (22) 

J (2 In/i 

Below, we will write eonvolution integrals of this form using the short-hand notation y = -{dZ/d In n)® 
Z“'. In the MS seheme, we have 


Z(m, 0 )',ijl) = 6{(jj - (o') + ^ ~ Z^^\co, (o',ju). (23) 

k=i ^ 

The Z faetors depend on fi both implieitly via a^f//), and explieitly via so-ealled star distributions 
related to the presenee of Sudakov double logarithms. The latter dependenee is a new feature, whieh 
leads to modifieations of the standard relations derived, e.g., in ina. Indeed, from the definition (El 
of the anomalous dimension it follows that 
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Here dasldlnju = /3(as) - 'leas is the generalized jS-funetion in the regularized theory, and P{as) is 
the ordinary yS-funetion. Both y and /3 are independent of e. Comparing ooeffieients of 1 /e^, we then 
find 

y = 2Qf, —— , (25) 
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While (El is a familiar result lEl, the relations (El eontain the additional dZ^''^ jd In /i pieee, whieh 
is usually not present. 

We now use the faet that, to all orders in perturbation theory, the anomalous dimension of the 
shape funetion is given by lfT!ni2l 
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where f{oS) is a smooth test funetion. It follows from (E51) that 
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cusp 


is the eusp anomalous dimension already mentioned in eonneetion with ©, whose two-loop 
expression has been derived in IfTKI . The other anomalous dimension, y, has been ealeulated at 
two-loop order in ll^ E3 • The relevant expansion eoeffieients are 
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The relations (EH) now allow us to express the eoeffieients in terms of the expansion eoeffi- 
eients of /3, Tcusp, and y. To derive these results, we need the following identities for star distributions 
(the first of whieh is somewhat laborious to derive) 
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Denoting by Z[„] the eoeffieient of (a,/An)" in Z(co, co',iu), we obtain after some algebra 
Z[o] = 6(co - co '), 
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This is the complete two-loop result for the renormalization factor of the 5-meson shape function. 

According to ©, the soft function is defined as the integral over the renormalized (parton-model) 
shape function. Using the fact that Z(m, only depends on the difference (co - co'), we find that 

Q, 

s{ — ,iu)= I JmZ(Q,m,yU)5bareM, (34) 
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where 5bare(f^) is defined as the integral over the bare shape function and is scale independent. Ex¬ 
panding this relation in perturbation theory, we obtain 
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1. The first term on the right-hand side in each line corresponds to the result obtained 
from the loop diagrams, given in (fTTll . The remaining terms correspond to operator counter-terms. 
Explicitly, we obtain for the counter-term contributions 
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where H^ 2 e is the harmonic number, which results form the integral 
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The counter-term contributions can be evaluated using the results for the bare one-loop soft func¬ 
tion from (El) and the expressions for the anomalous-dimension coefficients given in (EB- When 
adding these contributions to the result (El) for the bare soft function we find that all 1 /e" pole terms 
cancel, so that the limit 6^0 can now be taken. In O the logarithmic terms in the renormalized 
soft function have been determined by solving the renormalization-group equation for the function 
5. At two-loop order, it was found that 
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where the non-logarithmie one-loop eoefheient reads ifT^ BUI 



(39) 


Our results for the logarithmie terms agree with and thus eonfirm the existing results for the 
two-loop anomalous dimensions Fi IfHO and 71 lESlE^ll . In addition, our ealeulation gives for the 
non-logarithmie pieee at two-loop order the expression 


.( 2 ) 


Ct7 I--h 32^31 -l- CfC^ I- 

3 40 ^ / \ 81 

^ ^ ,8 5;r2 20 \ 


427;r2 

108 


61 n^ 



(40) 


This is the main result of the present work. 

It is interesting to eompare the exaet answer for the eoeffieient c® with the approximation ob¬ 
tained by keeping only the terms of order whieh eould be derived without any of the elaborate 

teehnology developed in the present paper. In the absenee of exaet two-loop results, it is sometimes 
argued that the /Soa^ terms eonstitute the dominant part of the eomplete two-loop eorreetion. In the 
present ease, we obtain for Nc = 3 eolors (note that /3o = 9 for = 3 light flavors) 

c® « 8.481 • j - 62.682 « -54.201. (41) 

Keeping only the term would give 8.481, whieh has the wrong sign and is off by almost an 
order of magnitude. This illustrates the importanee of performing exaet two-loop ealeulations. 


4 Discussion and summary 


Having eompleted the two-loop ealeulation of the soft funetion, we now briefly diseuss the impaet 
of our results for the predietion of the partial inelusive B X^y branehing ratio. We begin by 
displaying the final expressions for the funetions s{L, jx) and T(L, /i) obtained from (EH1> and © for 
the ease of 77c = 3 colors and n/ = 3 light quark flavors. We obtain 


s{L,n) ~ 1 + (-0.175 - 0.424L - OAIAL^) 

+ (-0.343 - 0.201L - 0.433L2 + 0.383L^ + 0.090L^) a]{pL) + ... , 

« 1 + (-0.873 - 0.424L - O.AIAL^) 

+ (-0.660 + 0.821L + 0.456L2 + + O.OOOL^) a]{pL) + ... . (42) 


The two-loop eorreetions are quite signifieant, espeeially in the ease of T(L,//). Figure |21 shows the 
dependenee of the soft funetions on L = ln(f 2 ///) at the fixed renormalization seale /i ehosen sueh 
that UsUi) = 0.45, eorresponding to a renormalization point n « 1.1 GeV as appropriate for the 
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Figure 2: One- and two-loop corrections to the soft functions s{L, jj) and T(L, /r) evaluated at a sip) = 
0.45. The dashed lines show the one-loop results results, while the solid lines give the complete two- 
loop results derived in the present work. The gray lines are obtained if only the terms are kept 
in the two-loop contributions. 


photon-energy cuts used in current measurements of the B Xsj decay. In addition to the one- 
and two-loop results, we display the results obtained if only terms of order are kept in the 
two-loop coefficients. The figure shows that the two-loop effects calculated in this paper can have 
an impact on the soft functions at the 10-20% level, and that keeping only the /3o«? terms does, 
in general, not provide an accurate description of the two-loop effects. We stress, however, that 
the large perturbative corrections seen in the figure do not translate in similarly large corrections 
to the B Xsj decay rate. The size of the corrections is strongly reduced once the pole mass in 
O is eliminated in favor of a low-scale subtracted h-quark mass, such as the shape-function mass 
ifT^ [T71 . At one-loop order this was demonstrated in IIT!n . and we expect similar cancelations to 
persist in higher orders. Note also that the soft functions by themselves are not renormalization- 
group invariant, so it is meaningless to study their dependence on the scale yu for fixed fl. In physical 
results such as the expression for the B Xsj decay rate in dU), the scale dependence of the soft 
function cancels against the /tq dependence of the objects U 2 , (A/yUo)'', and rj in ©. A detailed 
analysis of the phenomenological impact of NNLO corrections on the B Xgj decay rate will be 
given elsewhere. 

In summary, we have calculated the two-loop expression for soft function siL,fj.), which is de¬ 
fined in terms of an integral over the 5-meson shape function in the parton model. This quantity is a 
necessary ingredient for the NNLO evaluation of the B decay rate with a cut on the photon 

energy. For a sufficiently low cut energy, the partial inclusive decay rate can be calculated in a multi¬ 
scale operator-product expansion, in which the soft function arises in the final expansion step, when 
a current-current correlator in soft-collinear effective theory is matched onto bilocal heavy-quark 
operators in heavy-quark effective theory. If the cut on the photon energy is so severe that the con¬ 
tribution from the soft region cannot be evaluated perturbatively, the soft part should be subtracted 
from the partonic result for the cut rate, before it is convoluted with the renormalized shape function. 
Even in that case our result is a necessary component in the consistent calculation of the B X^y 
photon spectrum at NNLO. Moreover, since the soft function is universal to all inclusive heavy-to- 
light decays in the end-point region, our results are also relevant for semi-leptonic B XJ^v decay 
spectra, if one attempts to extend the analysis of to NNLO. 
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